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p^ Abstract 



A generalized kinetic model equation which takes into account the frequency depence of the thermal 
conductivity is used to analyze the problem of sound propagation in dilute polyatomic gases. By 
comparing the theoretical results with some available experimental data we infer that our model 
^N equation provides a precise transition between low and high-frequency limits. 
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^ 1 Introduction 

fj Essentially, there are two different approaches to solve the problem concerning the propagation of 

"r? plane harmonic waves in gases, namely: macroscopic and microscopic approaches. The macroscopic 

^;-> approach is based on conservation equations of mass, momentum and energy, and on the laws of 

Navier-Stokes and Fourier. Sound propagation predictions derived from this macroscopic approach 

are valid as long as the oscillation frequency of the acoustic wave is smaller than the mean molecular 

collision frequency [1]. In this case, sound dispersion is negligible, while the absorption per wavelength 

J> is proportional to the oscillation frequency and can be written as a sum of constributions due to viscous 

^^ and thermal effects. Different generalizations of the usual macroscopic approach - trying to extend 

^jl its validity to the high-frequency region - are reported in the literature [2, 3, 4]. Among them, the 

(Sj most popular are those that take into account the frequency dependence of the transport coefficients. 

.^ However, the sound propagation results derived from this approach give good agreement with the 

C^ experimental data only for the phase velocity of the sound waves. When the oscillation frequency is 

^^ comparable to the molecular collision frequency one has no recourse unless to turn to a microscopic 

• • approach, i.e., one has to use the Boltzmann equation. One of the earliest attempts to solve the 

. I^ problem of sound propagation in polyatomic gases based on a kinetic equation was made by Hanson, 

/\ Morse and Sirovich [5]. By employing the method of Sirovich and Thurber for polyatomic kinetic 

3 models - which is based on a generalization of the Gross-Jackson procedure for monatomic gases - 

they derived dispersion relations for the four-moment and seven-moment approximations. Calculations 

performed by Hanson, Morse and Sirovich for nitrogen and oxygen show that their results are in some 

agreement with the experimental results of Greenspan [6]. 

Our purpose in this paper is to study the propagation of sound waves in dilute polyatomic gases 
using a kinetic model equation which replaces the collision operator of the Boltzmann equation by a 
single relaxation time term [7] and generalizes the classical hydrodynamic description by considering 
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a frequency-dependent thermal conductivity. In particular, our expression for the generalized thermal 
conductivity follows from the Maxwell-Cattaneo law of heat conduction which can be written as a 
constitutive equation for the heat flux vector with an exponential memory kernel. By applying the 
normal mode method to our generalized model equation, it is possible to derive a dispersion relation 
which can be used to determine the phase velocity and the attenuation coefficient of the acosutic 
mode. Comparison of the theoretical results with the acoustic measurements performed by Greenspan 
in oxygen and nitrogen shows that our generalized kinetic model equation for dilute polyatomic gases 
provides a precise transition between low and high-frequency limits. 

We organize the paper as follows: in Section 2 we derive a generalized kinetic model equation for 
dilute polyatomic gases which takes into account the frequency dependence of the thermal conductivity. 
The sound wave propagation problem is analyzed in Section 3, while in Section 4 we compare the 
theoretical results with some available experimental data. Finally, in Section 5, we finish with some 
concluding remarks. 

Cartesian notation for tensors with the usual summation convention is used. Furthermore, angular 
parentheses around indices denote traceless symmetrization. 

2 Generalized Model Equation 

In the classical kinetic theory developed by Taxman [8] for dilute polyatomic gases, the one-particle 
distribution function f{x,c,t,s) is defined in such a way that f{x,c,t,s)dcds gives the number of 
molecules at position x and time t with molecular velocities between c and c + dc and with internal 
degrees of freedom between s and s + ds. By neglecting external forces, the one-particle distribution 
function satisfies the Boltzmann equation [9] 

where C{f, f) is the Boltzmann collision operator. The collision operator describes the rate of change 
of the distribution function due to molecular collisions and obeys the following collisional invariant 
conditions [10] : 



mC{fJ)dcds = 0, (2) 

l'mc^Cif,f)dcds = 0, (3) 



mc^ 



+ £s]C{f,f)dcds = 0, (4) 



where m is the molecular mass and Eg denotes the energy associated with the internal degrees of 
freedom. Conditions (2)-(4) correspond to the conservation of mass, momentum and energy during 
collisions, respectively. 

The mathematical complexity of the collision operator C(/, /) to analyze time-dependent problems 
like sound propagation and light-scattering in dilute polyatomic gases is usualy avoided by replacing 
it by a single relaxation-time term of the form 

C{fJ) = -^-^^^, (5) 



where r is an effective relaxation time and fr is a reference distribution function which satisfies the 
main physical properties of the Boltzmann collision operator. Recently, an expression for the reference 
distribution function was obtained by Marques Jr. [11] by requiring the Chapman-Enskog solution of 
the kinetic model equation to be consistent with the classical Navier-Stokes-Fourier description. In 
the usual Navier-Stokes-Fourier theory a macroscopic state of a dilute polyatomic gas is characterized 
by the fields of mass density 



flow velocity 



mfdcds, 



(6) 



mcifdcds, 



and temperature 



while the pressure tensor 



and the heat flux vector 



ipcvr' I ( 



mC^ 



+ £s I fdcds, 



Pij 



I 



mCiCj fdcds 



mC^ 



+ £s ) Ci fdcds 



are given, respectively, by the following constitutive relations [12] 



Pij 



dVr 



dv 



^-^»9^j^^^-'^"ax 



{i 



'i> 



and 
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-Ao 



dT 



(7) 



(8) 



(9) 



(10) 



(11) 



fl2) 



where c^, is the total specific heat at constant volume, Ci = Ci — vi is the peculiar velocity, p is the gas 
pressure, ^o is the shear viscosity, r]o is the bulk viscosity and Aq is the thermal conductivity. By taking 
the effective relaxation time r equal to the stress relation time r^ = fJ^o/p, the reference distribution 
function for dilute polyatomic gases which is compatible with the classical Navier-Stokes-Fourier theory 
reads 



/, = /(") 1 + rl 



fo 



mC^ 



£ 



^U(^ 

2kT 2 \kT kT 



adT 

T dxi 



+ t(1-Zo) 



where 
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kT kT 



^(0) _ r_ 



p/ m \3/2 (mC^\e-x:p{-£s/kT) 



dVr 
dXr 



(13) 



(14) 



is the equilibrium distribution function and 

Z= I e^p{-£s/kT)ds (15) 

is the partition function associated with the internal degrees of freedom. Moreover, 7 is the specific 
heat ratio, f^ = Ao//U(,Ct, is the Eucken factor, Z^^ = (ryo//Uo)/(5/3 — 7) is the internal relaxation number 
(i.e., the mean number of molecular collisions required for the translational and internal degrees of 
freedom to come to thermal equilibrium [13]) and 

£=^ fSsexp {-£s/kT) ds (16) 

is the mean internal energy. 

By analyzing time-dependent problems like sound propagation and light scattering in dilute poly- 
atomic gases. Marques Jr. [11] was able to determine the range of applicability of his kinetic model 
equation. Comparison of theoretical results with available experimental data in nitrogen, oxygen, 
carbon dioxide and methane shows that the reference distribution function (13) is valid as long as the 
external oscillation frequency is smaller than the relaxation frequency for internal and translational 
degrees of freedom to come to thermal equilibrium. 

It is possible to extend the validity of the kinetic model equation to the high-frequency region by 
considering a generalization of the classical hydrodynamic description. Certainly, the most popular 
approach that generalizes the hydrodynamic behaviour of simple fluids is based on the Maxwell- 
Cattaneo law of heat conduction [3] 

where the relaxation time r, gives us a measure of the time interval spent by the heat flux vector to 
achieve a stationary value. It is well known [14, 15] that the Maxwell- Cattaneo equation (17) can be 
written as 

qi = - -expl — o . dt, (18) 

JO ^1 \ ^q J O'Xi 

i.e., as a constitutive equation with an exponential memory kernel in such a way that a generalized 
thermal conductivity can be defined as 

A(t-fO = -expf-^^V (19) 



Since sound wave and light scattering solutions are found by Fourier transformation in space and time, 
we verify from (19) that this generalization introduces a frequency-dependent thermal conductivity 

A(w) = r-^ . (20) 

Furthermore, based on the work of Zwanzig [4], it is possible to modify the constitutive relation 
(11) for the pressure tensor to take into account the frequency dependence of the bulk viscosity. 
Aiming to calculate the spectral distibution of scattered light in a one-component fiuid whose molecules 
have internal degrees of freedom weakly coupled to their translational degrees of freedom by a single 
relaxation time process. Mountain wrote the pressure tensor as [16] 



where r* is the relaxation time. Expression (21) shows that in the so-cahed weak couphng hmit the 
bulk viscosity consists of two parts, a frequency independent one due to translational motions and a 
frequency dependent one which is related to the exponential decay of the energy in internal modes as 
a result of interactions with translational modes. Comparison between theory and experiments shows 
that the generalized hydrodynamical model proposed by Mountain does not always apply since in 
some fluids the relaxation of the internal degrees of freedom involves more than one relaxation time. 
Based on these facts we construct in this work a generalized kinetic model equation for dilute 
polyatomic gases which takes into account the frequency dependence of the thermal conductivity via 
the Maxwell-Cattaneo law, but disregards the frequency dependent part of the bulk viscosity. We 
start by assuming that the reference distribution function is given by the expression 



fr = f 
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where A, Ai, Arr, ^(ij) and Airr are expansion coefficients that depend on position and time through 
the basic fields. As described in details by Marques Jr. [11], the application of the Chapman-Enskog 
method [9] to solve the kinetic model equation allows the determination of the expansion coefficients 
appearing in the above reference distribution function. Hence, by considering the constitutive relations 
(11) and (18) we get 
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Finally, by combining expressions (1), (5), (22) and (23) we can write our generalized kinetic model 
equation for dilute polyatomic gases as 

'mC^ 3\ , 
7 ) I ITTT^ - o - (7 



dt * dxi 
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is the generalized Eucken factor. 

Closing this section, we call attention to the fact that the application of the above generalized 
kinetic model equation to study time-dependent problems in dilute polyatomic gases only requires 
the specification of the ratio of the specific heats 7, the Eucken factor fg and the internal relaxation 
number Zq. 



3 Absorption and Dispersion of Sound 

As an application of the generalized kinetic model equation derived in the previous section, we study 
the problem concerning the propagation of a plane harmonic wave in dilute polyatomic gases. To 
simplify this study we consider the sound propagation process in the linear regime near equilibrium 
and assume that the sound wave moves along the x-axis. Hence, we write the distribution function as 

f = fo{l + i^exp[i{Kx-ojt)]}, (26) 

where /o is the absolute equilibrium distribution function, k = oj/v + ia is the complex wavenumber, 
CO is angular oscillation frequency, v is the phase velocity, a is the attenuation coefficient and V' = tp{c) 
is a function that depends only on the molecular velocity. In sound propagation problems the complex 
wavenumber k is determined as a function of the oscillation frequency uj via solution of the dispersion 
relation (for details see [11, 17, 18]) 

det (7^M(°^ - i (kvo/uj) (I - M'^^)) = 0, (27) 

where 7^ = po/fio^ is the rarefaction parameter, v^ = y^2kT^/m is the equilibrium thermal velocity 
and I is the identity matrix. Moreover, M*''' and M'^' are 3x3 matrices whose elements - given in the 
Appendix - depend on the dimensioless parameter z = {u}/kVu){1 + iTZ) and on the plasma dispersion 
function [19] 

1 /•+°^ exp (-t^) 
W(z) = -^ / -^-^ ^-dt. 28 



The acoustic solution of the dispersion relation (27) in the low-frequency limit can be determined 
by expanding the dimensionless wave number kvq/oj in powers of 1/7^. By retaining terms up to first 
order we obtain 



KCq 

UJ 


-'+27W 
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(29) 



We verify from expression (29) that the dispersion of sound is negligible in the low-frequency region, 
while the sound absorption per wavelength is proportional to the sound oscillation frequency and can 
be written as a sum of viscous and thermal effects. Note that an additional contribution coming from 
internal relaxation appears in the absorption coefficient of polyatomic gases, a fact that allows us to 
determine the bulk viscosity by comparing the theoretical results with available experimental data. 

4 Comparison with Experiments 

In order to compare the sound propagation results derived from our generalized kinetic model equa- 
tion (24) with the acoustic measurements performed by Greenspan [6] on nitrogen and oxygen we shall 
use the following material parameters: 

Gas 7 fo Zq 

^^^ L40 L97 2^ 

O2 1.39 1.95 1.3 



The values of the specific heat ratio 7 and the Eucken factor fg were obtained by using the values 
of c„, Cp, ^0 and Aq given in the CRC Handbook of Chemistry and Physics [20], while the values of 
the internal relaxation number Zq were obtained by using expression (29) to fit the sound absorption 
experiments of Greenspan in the low- frequency limit. 

In Fig. 1 and Fig. 2 the attenuation factor aco/co and the reciprocal speed ratio Cq/v are shown 
on a double logarithmic scale as a function of the rarefaction parameter TZ = po/fiooo for nitrogen and 
oxygen, respectively. The solid line represents the theoretical sound propagation results derived from 
our generalized kinetic model equation, while the dashed line represents the theoretical results derived 
from the usual macroscopic (hydrodynamic) approach based on the laws of Navier-Stokes and Fourier. 
The open circles are the experimental data of Greenspan for sound wave absorption (lower curve) and 
dispersion (upper curve). 
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Figure 1: Attenuation factor aco/uj and reciprocal speed ratio Cq/v as a function of the rarefac- 
tion parameter TZ = p^/ ^^u for nitrogen at 300K. The predictions of our generalized kinetic model 
equation ( — ) are compared with the experimental data of Greenspan and the results of the usual 
macroscopic approach ( — ). 



We verify from Figs. 1 and 2 that in the low-frequeny region [TZ ^ 1) the theoretical sound 
propagation results derived from our generalized kinetic model equation are in complete agreement 
with the experimental data. Moreover, the predictions for sound wave absorption and dispersion are 
exactly the same as obtained by the macroscopic approach. Since we required the Chapman-Enskog 
solution of our model equation to be consistent with the usual macroscopic description, this result is 
in complete agreement with our expectations. In the transition region (1 < 7^ < 10) we note that 
the kinetic theory proposed in this paper yields better results than the usual macroscopic approach. 
Whereas the latter predicts a value for the attenuation factor which is about 23% greater than the 
experimental value, the predictions derived from our generalized kinetic model equation are in good 
agreement with the absorption data up to a value of the rarefaction parameter of approximately 2. 
For smaller values of the rarefaction parameter, i.e., in the high-frequency region we observe that 
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Figure 2: Attenuation factor aco/uj and reciprocal speed ratio Cq/v as a function of the rarefac- 
tion parameter TZ = P(,/fioUj for oxygen at 300 K. The predictions of our generahzed kinetic model 
equation ( — ) are compared with the experimental data of Greenspan and the results of the usual 
macroscopic approach ( — ). 



the usual macroscopic approach yields a qualitatively correct result for the absorption curve but fails 
to describe the dispersion data. In contrast, the theoretical results obtained from our generalized 
kinetic model equation for small TZ are in good agreement with both, especially with the sound wave 
dispersion curve of nitrogen. Further, we can conclude from Figs. 1 and 2 that the agreement between 
theory and experiment in the high-frequency limit is better for nitrogen than for oxygen, i.e., for gases 
with higher values of Zq which can be explained in terms of the internal relaxation time r = HoZo/p. 
A higher value of Zq means that the thermal relaxation process between translational and internal 
degrees of freedom runs more slowly and, therefore, has a greater influence on the propagating sound 
wave which is appropriately described by our generalized kinetic model equation. 

The comparison presented in Figs. 1 and 2 indicates that there are still opportunities for fur- 
ther enhancements in the kinetic theory, such as including the complete frequency dependence of all 
transport coefficients or considering non-linearities in the model equation, but all in all, we want to 
remark that the kinetic theory is well suited to describe sound propagation in dilute polyatomic gases 
over a wide range of frequencies and provides a precise transition between low and high-frequency 
limits. Whereas the macroscopic approach becomes unable to describe the acoustic measurements in 
the high-frequency region, the kinetic model equation derived in this work still provides a qualitatively 
correct description of the experimental data, even for values of the rarefaction parameter smaller than 
the unity. 



5 Conclusion and Outlook 

In the present work, we are concerned with the time-dependent problem of sound propagation in dilute 
polyatomic gases with classical internal degrees of freedom. Since the macroscopic (hydrodynamic) 
approach fails to describe the experimental data in the high-frequency region one needs to turn to 
a microscopic (kinetic) approach based on the Boltzmann equation. Assuming that the Boltzmann 
collision operator can be replaced by a single relaxation-time term in order to eliminate its mathemat- 
ical complexity, we are able to derive a generalized kinetic model equation for dilute polyatomic gases 
which explicitly takes into account the frequency dependence of the thermal conductivity. By com- 
paring the theoretical sound propagation results derived from our generalized kinetic model equation 
with the experimental data, we observe that the kinetic theory is well suited to describe sound propa- 
gation in dilute polyatomic gases over a wide range of frequencies. Whereas the macroscopic approach 
becomes unable to describe the acoustic measurements in the high-frequency region, the generalized 
kinetic model equation proposed in this work still provides a qualitatively correct description of the 
experimental data, even for values of the rarefaction parameter smaller than the unity. The presented 
comparison also indicates that there are still opportunities for further enhancements in the kinetic 
theory. One of the next steps for future work on this topic would be to include the complete frequency 
dependence of all transport coefficients into the kinetic model equation. 

Finally, we want to close this section by remarking that for the application of our generalized 
kinetic model equation to a specific sound propagation problem no further information about the 
potential energy of interaction between the gas molecules is required because the acoustic properties 
are completely characterized by three parameters, namely the ratio of the specific heats 7, the Eucken 
factor i(, and the internal relaxation number Zq. Values of these parameters can easily be found in the 
literature or fitted to experimental data in the low- frequency limit. 



Appendix 

The non-zero elements of the matrices M'''^ and M^^^ read 



MS' 



M' 



(0) 



) ^•■^22 



M(°'=zM<«', M<«) = zMg) 
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(31) 



M(«) = (7 - 1)M(°3), MW = 2(7 - 1)M<«3), Mg) = (7 - l)D{z) + ^ (^ " t) W{z), 
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(36) 



being the frequency-dependent Eucken factor. Moreover, we have introduced the abbreviations 

A{z) = zW{z) + l, B{z)=(z^-]^W{z) + z, C{z) = z(z'^-^W{z) + z^-l, (37) 



D{z) =[z'^-z'^ + -\ W{z) + z^ 



% E{z)=U-z^-^^W{z)^z^-\, (35 



F{z) = z{z^-2z^ + -] W{z) + z^--z'' + -. 



(39) 
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